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SUMMARY 


Fracture mechanics has been found to be a useful tool for understanding 
composite delamination# Analyses for calculating strain energy release rates 
associated with delamination growth have been developed# These analyses 
successfully characterized delamination onset and growth for particular 
sources of de lamination# Low velocity impact has been found to be the most 
severe source of composite delamination. A variety of test methods for mea- 
suring interlaminar fracture toughness are being developed to identify new 
composite materials with enhanced delamination resistance. 

INTRODUCTION 

As high strength, low density fiber reinforced polymer matrix composites 
are considered for highly strained primary aircraft structures, increased 
attention is being devoted to the understanding and characterization of sig- 
nificant composite failure modes# Perhaps the most commonly observed failure 
mode in composites is delamination, which is a separation of the individual 
plies# Delamination may result from low velocity impact, from eccentricities 
in structural load paths that induce out-of -plane loads, or from discontinu- 
ities in the structure which create local out-of-plane loads# Because delam- 
ination occurs so frequently, interest in analyzing and characterizing delam- 
ination behavior has increased# Fracture mechanics has been found to be a 
useful tool for this purpose# Summaries of some o* the significant recent 
studies on delamination are included in this review# 

A Perspective 

There is a fundamental difference between the way linear elastic fracture 
mechanics is being applied to delamination in composites and the way it is 
used for metals. In metals, the largest crack in a component is assumed to 
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control strength* Failure of the structure may be directly related, by the 
material^ fracture toughness, to catastrophic propagation of this crack. 
However, composite delamination in structures subjected to in-plane loading is 
a subcritical failure mode whose effect may be (1) a stiffness loss that is 
benign in terms of the structural failure, (2) a local tensile strain concen- 
tration in the load bearing plies that causes tensile failure, or (3) a local 
instability that causes further growth which ends in compressive failure. In 
the latter two cases the delamination leads to a redistribution of structural 
load paths which, in turn, precipitates structural failure. Hence, de lamina- 
tion is indirectly responsible for the final failure of a composite, whereas 
crack propagation in a metal is directly responsible for structural failure. 
Consequently, interlaminar fracture toughness for composite delamination does 
not have the same significance to design as fracture toughness does for a 
metal. The designer must be able to determine the consequence of the delami- 
nation in terms of structural load redistribution and relate this to some 
appropriate structural failure criterion. Nevertheless, fracture mechanics is 
a useful tool for understanding the mechanics of delamination, to determine 
the variables that control when delaminations form and grow, and for charac- 
terizing the inherent delamination resistance of the composite. This paper 
will focus on these aspect? of the fracture mechanics characterization of 
delamination in composite materials. 

MECHANICS OF DELAMINATION 

Strain Energy Release Rate Analysis 

Some potential delamination sites are the common design details that 
result in discontinuities in the load path. These discontinuities give rise 
to interlaminar stresses even when the remote loading is in-plane [1]. Some 
of these typical design details (fig* 1) are (1) a straight edge with a 
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stress-free boundary (often called a "free edge"), (2) a curved edge such as a 
hole boundary, (3) a drop off of the interior plies of a laminate to taper 
thickness, (4) a bonded or co-cured joint, and (5) a bolted joint. Of these 
cases, the free-edge delamination problem has been studied most extensively, 
and illustrates the benefits of fracture mechanics for characterizing 
delamination behavior. 

Several investigators have analyzed the strain energy release rate, G, 

associated with edge delamination growth [2,3,4]. From finite element analysis 

they have found that once the delamination progresses beyond a distance equal 

to a few ply thicknesses from the edge, G reaches a constant plateau given by 

the first equation in figure 2. This equation, derived from a rule of mixtures 

and laminated plate theory, shows that G is independent of delamination size 

and depends only on the remote strain, e, the laminate thickness, t, and two 

modulus terms, E LAM an< * E *' that correspond to the laminate modulus before 

and after delamination. In Reference [4], the critical strain, e c , measured 

at the onset of edge delamination in the -30/90 interfaces of eleven-ply 

[± 30/± 30/90/90] laminates was substituted into the edge delamination equacion 
s 

in figure 2 to determine a critical G c for Thornel 300/Narraco 5208* 

(graphite/epoxy) composites. This G c was then used to predict onset strains 

for edge delamination in the 0/90 interfaces of [+45 /-45 /0 /90 ] laminates 

n n n ns 

of the same material, where n * 1,2,3 corresponded to and 8-, 16-, and 24-ply 

laminate. Figure 3 shows that the predictions agreed well with experimental 

data, and captuied the trend of decreasing delamination onset strain with 

increasing thickness. An analysis based on critical values of interlaminar 

stresses would not account for this thickness dependence because the 

interlaminar stresses are independent of n [3]. Ho rigorous analytical 

explanation has been presented to explain the ability of the strain energy 
*The use of trade names in this article does not consitute endorsement, 
either expressed or implied, by the National Aeronautics and Space 
Administration or AVSCOM. * 
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release rate (which is rigorously defined only when a crack of finite size 
exists) to predict delamination onset when a stress criterion cannot. 
Although the interlaminar stresses in the interfaces that delaminate are 
singular at the straight edge, these si; gularities are not the same as 
those that exist when the delamination has formed [5,6]. Perhaps the most 
logical explanation is the presumption that inherent interlaminar flaws of 
microscopic dimension are initially present in the interface near the edge, 
or they are created by the singular stresses at the edge, but they do not 
form delaminations of a finite detectable size until a critical G value 
is reached for delamination extension [4,7,8]. 

A similar analysis for a local de lamination growing from a matrix ply 
crack in an off-axis ply yielded the second equation shown in figure 2 
[9]. This equation also shows that G is independent of delamination size 
and depends only on the remote load, P, the laminate width, W, the lami- 
nate thickness, t, the laminate modulus, and the thickness and 

modulus of the uncracked, locally delaminated regions, t^j) and E LD , 
respectively. Pigure 4 shows predictions (open symbols) of delamination 
onset in the -25/90 interfaces of [±25/90 n ] s laminates. Under tensile 
loading, these laminates undergo a transition from edge delamination 
(n < 3) to local delamination growing from matrix cracks (n > 4 ). A criti- 
cal G c was calculated from the n * 2 edge-delamination data, and then 
was substituted into the appropriate equations to predict delamination on- 
set strains for the other layups. As shown in figure 4 , the G anlyses 
predict the onset of both types of delamination well. Furthermore, just as 
the dependence for edge delamination was predicted by G analysis (fig. 

3), so the trend of decreasing e c with increasing thickness of the 
cracked ply, and hence increasing n, was predicted for local delamination 
(fig. 4). 
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The ability of the strain energy release rate to correlate delamination 
behavior from several sources and account for volumetric effects in the form 
of thickness dependence is the primary reason for adopting fracture mechanics 
to analyze delamination problems. Historically, this same capability for 
predicting crack growth in large metallic structures based on fracture tough- 
ness measurements of small test specimens led to the acceptance and growth of 
fracture mechanics in the 1950's [10]. 

Mixed-Mode Behavior 

Because a delamination is constrained to grow between individual plies, 
both interlaminar tension and shear stresses are commonly present at the 
delamination front. Therefore, delamination is often a mixed-mode fracture 
process. A boundary value problem must be formulated and solved to determine 
the Gj f Gjj, and Gju components of the total strain energy release rate. 

A virtual crack extension technique has been combined with displacement- based 
finite element analysis to calculate the various G components [2]. This 
technique, in conjunction with quasi-three-dimensional finite element anal- 
yses, was applied to the free-edge delamination problem by several authors 

[2.3.4.7.11.12.13] . Figure 5 shows the results of such an analysis for the 
[ ±30/±30/90/90 ] laminate when each of the -30/90 interfaces are assumed to 
delaminate. As noted earlier, the total G raises to a constant plateau when 
the delamination progresses beyond a few ply thicknesses from the edge. 

Figure 5 indicates that the total G is equal to the sum of Gj and Gjj 
components. The Gjjj component calculated was negligible. Furthermore, 
like the total G, the Gj and G IZ components are also independent of 
delamination size for edge delamination. 

Similar analyses have been performed for other composite layups 

[3.11.12.13] . Figure 6 shows the Gg percentage of the total G for three 
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quasi-isotropic layups containing ±45°, 0 3 , and 90° plies. Each layup had 90° 
plies in the center, bit the order of the +45°, -45°, and 0° plies was permu- 
tated. De laminations were modeled between the 90° plies and adjacent angle 
plies. The total G was nearly identical for the three layups but, as shown 
in figure 6, the Gj component varied significantly. 

Because of the complex mixed-mode nature of composite delamination, no 
closed form solutions have been developed to lay a strong theoretical founda- 
tion for understanding the governing parameters that control delamination 
behavior. Recently, however, a singular integral equation approach based on 
anisotropic composite laminate elasticity theory was used to establish some 
general asymptotic solutions for edge delamination [14]. A partially closed 
crack model for delamination was introduced to remove the restriction of 
physically inadmissible oscillatory stress singularities. Both the global 
response in terms of strain energy release rates (Gj, G jl' g ih' an< * 

G*iOTAL ) and tiie asymptotic field solution in terms of crack-tip stress inten- 
sity factors (Kj, K XI , and Kjjj) were calculated. Anisotropic laminate 
elasticity solutions were obtained for relatively simple cases of delamination 

in (+0/-0) laminates [14]. Finite element methods based on a combined 

s 

laminate elasticity solution and singular hybrid finite element formulation 
were employed for fiber composites with complex lamination variables [15]. 
Figure 6 shows the Gj percentage of the total G for edge delamination in 
the ( ±45/0/90 ) g family analyzed previously. The same trend in Gj/G with 
layup is apparent, indicating that the simple approximate techniques, used 
with the more common displacement-based finite elements, may be reasonably 
accurate for predicting delamination behavior. The virtual crack extension 
method has been successfully used with a 3D NASTRAN finite element analysis to 
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predict delamination growth from both a manufacturing flaw and an open hole in 
the F-16 composite horizontal tail [l]. 

The fracture mechanics analyses of other sources of delamination will 
continue to be a fruitful research topic as composite materials become 
candidates for primary aircraft structure. Both the U.S. Air Force and the 
U «S* Army are developing damage tolerance and durability criteria for future 
composite structures which may ultimately require such analyses for structural 
certification* 

DELAMINATION GROWTH 

Delaminations in composite laminates may interact with other damage 
mechanisms that result in different growth behavior than may be anticipated 
from an elastic analysis* For example, under quasi -static loading, edge 
de lamination formation is followed by matrix cracking in the 90° plies* These 
matrix cracks may extend throughout the laminate width, ahead of the delamina- 
tion front [4], causing a perturbation in the three-dimensional stress field 
at the delamination front* This interaction between matrix cracks and delami- 
nation causes stable de lamination growth. A two-dimensional linear elastic 
fracture mechanics analysis cannot model this three-dimensional damage, but it 
may be used for semi-empirical characterization of stable growth using the 
R-curve method [16]* Figure 7 shows an R-curve for stable edge delamination 
growth in a [±45/0/90] o T300/5208 laminate [17]* Unlike the G predicted 
from an elastic analysis that is independent of delamination size, a, the 
delamination resistance, Gg, increases with a as the matrix cracking 
accumulates ahead of the delamination front. As shown in figure 8, this 
R-curve was used to predict stable edge-de lamination growth in a thicker 
l+452/-452/02/90 2 l s laminate, and in a t±30/±30/90/90] laminate with a dif- 
ferent layup* 
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If the composite has a tough matrix, however, the matrix cracking may 
be suppressed. Delaminations will form at higher G values, but they will 
immediately extend through the laminate width in an unstable manner (7]. 
This edge delamination behavior in a toughened-matrix composite is the 
opposite of the stable crack growth observed in a ductile metal where large 
plastic zones at the crack tip result in more stable growth. 

Another source of delamination resistance that results in stable 
growth is the fiber bridging that occurs when a delamination is forced to 
grow between plies of similar orientation. This often occurs in the double 
cantilever beam (DCB) test used to measure Mode I interlaminar fracture 
toughness [18,19,20]. Figure 9 shows the DCB specimen, which typically 
consists of 24 plies of unidirectional material. As shown on figure 9, the 
specimen compliance is a non-linear function of delamination size. In ref- 
erence [20], a Teflon* coated glass fabric was placed at the midplane in 
the end of the beam to simulate a crack-like sharp notch. Measurements of 
Gj c were taken at the initial insert length (delamination length ■ 0) and 
along the length of the beam as the delamination grew. Figure 10 shows the 
R-curve generated from a DCB specimen of Hercules HMS/3501-6* (graphite/ 
epoxy) . The three-fold increase in Gj c from the initiation at the crack 
starter until values stabilized after the delamination grew over one third 
of the beam legnth was attributed to fibers debonding from the matrix, 
bridging the delamination plane, flexing, and breaking. A similar R-curve 
generated for an Hercules AS-1/3501-6* DCB specimen showed a higher G Ic 
at initiation and less increase in G Ic as the delamination grew (fig. 

10). In this case, much less fiber bridging was observed due to the 
stronger interfacial bond between the AS-1 fiber and the 3501-6 matrix. 

Unlike the stable delamination growth that occurs under quasi-static 
loading at G levels above G c , cyclic delamination growth may occur at 
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levels well below the static interlaminar fracture toughness. Figure 11 shows 
results for cyclic edge delamination tests at maximum cyclic strain levels 
below the delamination onset strain measured in the static test [11]. Delami- 
nations formed at these lower cyclic strains after a certain number of cycles, 

N. G c values, calculated from maximum cyclic strains and plotted as a func- 
tion of the number of cycles to delamination onset, are shown in figure 11. 

These G c values drop sharply in the first 200,000 cycles until they reach a 
plateau tantamount to a threshold for delamination onset in fatigue. An 
elastic analysis of the [±35/0/90]_ and [0/±35/90] laminates yields an iden- 
tical total G for a given nominal strain, but their Gj components are very 
high and low, respectively. As shown in figure 11, the G Q for these two 
layups measured in static tests (N * 0) are different. However, the G c 
thresholds under cyclic loading are nearly identical. Hence, the total G 
does not determine the delamination onset under static loading, but does 
appear to govern the de lamination onset in fatigue. Furthermore, comparison 
of the relatively tough Hexcel HX205 matrix composites to the brittle 5208 matrix 

composite shows a significant improvement in the static G c , yet the magnitude 
of this improvement for the G c threshold in fatigue is much less. 

The comparisons shown in figure 11 are useful for selection of more 
durable composite matrices. However, for damage tolerance analysis, a com- 
plete characterization of elamination growth, from the onset threshold under 
cyclic loading, G t ^, to the static interlaminar fracture toughness, G c , is 
required (fig. 12). Several investigators have used the mode I DCB and mixed- 
mode cracked- lap-shear (CLS) specimens to characterize this delamination 
growth [1,21,22,23]. They have developed delamination growth laws of the form 
~ m cG b similar to cyclic crack growth laws used to characterize fatigue 
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crack growth in metals. These baseline data were used with finite element 
analyses to predict cyclic delamination growth in structural components [1]. 

One source of composite delamination growth is the local instability that 
arises from the buckling of a delaminated region under compression loading 
[24,25]. Figure 13 shows a comparison of measured delamination cyclic growth 
rates, da/dN, and strain energy release rate components, Gj and j , 
calculated from a geometrically nonlinear finite element analysis for a 
compressively loaded laminate containing a through-width delamination [25]. 
Both Gj and the delamination growth rate increase and then decrease rapidly 
with delamination extension. However, G^j increases monotonically to a 
constant level with delamination extension. Hence, it appears that Gj is 
primarily responsible for driving the de lamination in the buckled region, even 
though the Gjj component was an order of magnitude larger. 

Figure 13 also shows the tendency for the delamination cyclic growth rate 
to decrease with delamination length. Hence, the growth of a single delamina- 
tion under cyclic compression may arrest and hence may not severely degrade 
performance. Figure 14 shows plots of cycles to failure as a function of 
maximum cyclic strain for a 38-ply T300/5208 (graphite /epoxy) laminate, typi- 
cal of an upper surface transport wing skin, subjected to cyclic compressive 
loading [26] where (1) the laminates contain a single circular insert, of 
various diameters, located at different interfaces to simulate delamination, 
(2) the laminates have been subjected to low velocity impacts over a range of 
impact energies, and (3) the laminate contains a 1.27 cm (0.5 in) diameter 
hole. The impacted laminates suffered severe degradation in the cosipressive 
strain at failure. In most cases the static compressive failure strain after 
impact was less than the residual compressive strain for the implanted de lami- 
nation after 10® fatigue cycles. 
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Figure 15 shows a plot of cycles to failure as a function of stress 
amplitude for (0/90/0/±45/0]_ graphite epoxy laminates subjected to fully 
reversed cyclic loading^ either in the undamaged state^ or following an impact 
with a potential energy per unit thickness of 1790 J/m [27). All failures 
occurred in the compression cycle. The data in figure 15 indicate that the 
compressive strength after impact is very low compared to the fatigue behavior 
of the virgin laminate. Furthermore, most of the strength reduction occurs 
after the impact, with very little degradation due to subsequent cyclic 
loading. 

The results in figures 14 and 15 indicate th<_, unlike metals where 
cyclic loading poses the most severe threat of crack growth leading to fail- 
ure, the most severe source of delamination in graphite epoxy structural 
laminates is low velocity impact. Such impacts may create significant delami- 
nation throughout the thickness of the laminate, degrading compressive 
strength. Furthermore, this extensive delamination may occur on the laminate 
interior with little or no visible damage apparent on the laminate surface 
[26-30]. Some research has just begun to relate the strain energy release 
rate for delamination growth due to out-of-pla te loading during impact to the 
growth of multiple de laminations that form throughout the laminate thickness. 
One analysis is based on large deflection plate theovy for circular isotropic 
plates containing multiple delaminations [31]. The influence of plate thick- 
ness, support conditions, and matrix toughness on the onset and propagation of 
delamination was predicted. Figure 16 shows a comparison of the predicted and 
measured load-displacement record for a 16-ply, circular quasi-isotropic 
plate, 25 am in diameter, subjected to a transverse load. Seven symmetrically 
loaded delaminations were assumed to grow simultaneously. Good agreement was 
found for the initial delamination onset at a load of 1 ks, and for the 
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load-displacement curve, up to the point that fiber breakage occurred in 
the back surface ply, at a load of 3 kN. More refined analyses that model 
the specific damage details observed during low velocity impact are needed 
to further clarify the relationships between impact loading and interlami- 
nar fracture. 

CHARACTERIZATION OP DELAMINATION RESISTANCE 

Because of this severe delamination threat posed by non-visible, low 
velocity impact, there has been increased interest in developing new tough- 
ened matrix resins for graphite reinforced com^^sites. Several investiga- 
tors have identified a correlation between improved compressive strength 
after impact and improved interlaminar fracture toughness for oompositses 
with toughened matrix resins (26,29,30]. These observations provided moti- 
vation for developing new test methods for measuring interlaminar fracture 
toughness. 

Many investigators have used a double cantilever beam (DCB) specimen 
(fig. 9) to measure critical values of the crack opening (interlaminar 
tension) fracture toughness, Gf C [18-23,32-34]. The edge delamination 
test (EOT) (fig. 17) is also bein^ used to measure the mixed-mode fracture 
toughness of composites [7] and determine fracture mode dependence (35,- 
36]. Figure 18 shows a comparison of the interlaminar fracture toughness 
of graphite composites with very brittle or tough matrix resins tested 
using a 24-ply unidirectional DCB as well as three EOT layups with low, in- 
termediate, and high percentages of Gj . The solid symbols represent mode 
I component of measured G c . Results indicated that only the Gj compon- 
ent controlled interlaminar fracture for the brittle matrix composites 
(5208, 2202-1)* whereas both Interlaminar tension and shear contributed to 
delamination in tough matrix composites (HX205/Fl«*5f *. Both the DCB and 
EOT testa are being used to evaluate new composite materials developed for 
heavily loaded large composite wing skins under the RASA Aircraft Energy 
Efficiency (ACKE) project (37). 
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Several other tests are also being used to measure interlaminar fracture 
toughness, A pure mode I version of the EET test has been proposed [38], The 
cracked-lap shear-(CLS) test for mixed mode G c has been used by several 
investigators [19,21,22], Currently an ASTM task group is conducting a round 
robin test program using the DCB, EDT, and CLS specimen as a prelude to devel- 
oping ASTM standards. Other investigators have used unidirectional beams 

loaded in three-point bending to estimate g [19,39,40,41] or loaded 

lie 

asymmetrically to measure mixed mode G c values with different ratios of 
G l/ G H [19,40], 

The recognition of the importance of interlaminar fracture toughness to 
composite damage tolerance has led to the development of new toughened matrix 
composites. Because resin manufacturers would prefer to measure the fracture 
toughness of a small neat resin specimen, there has been renewed interest in 
composite microraechanics analyses. The ultimate goal of such analyses would 
be to relate the fracture behavior of the neat resin to the interlaminar frac- 
ture behavior of the composite. Such analyses would guide the development of 
neat resin fracture tests that would yield realistic estimates of improvements 
in interlaminar fracture toughness of composites, 

CONCLUSIONS AND FUTURE DIRECTIONS 

Although the promise of composite materials has been known for several 
decades, the application of these materials to primary aircraft structure is 
just beginning. As suitable toughened-resin composites are developed to 
overcome the existing impact problems, methods for characterizing their 
behavior will have to be re-evaluated. For example, it may become necessary 
to use a J-integral analysis to correlate delamination behavior between 
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different layups and structural discontinuities. Furthermore, fatigue may 
reappear as the limiting factor for damage tolerant design. In any event, the 
ability to anticipate and quantify delamination failure using fracture mechan- 
ics analyses will still play an important role in the understanding of 
composite behavior, in establishing the damage tolerance and durability design 
criteria for certification, and in the development of improved composite 


materials 
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Fig 2 Sources of delamination in unnotched laminates 




f +45/ -45 _/ 0 / 90 I GRAPHITE EPOXY 
L n n n nj$ 


ORIGINAL PAGE I# 
OF POOR QUALITY 





' 


NUMBER OF PLIES 

Fig 3 Edge delamination onset in 0/90 interfaces of T30Q/5208 laminates 
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Fig 5 Strain energy release rate components for edge delamination 
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Fig 6 Mode I percentage of total G for edge delamination in 6/90 interfaces of 
quasi -isotropic laminates 
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Fig 8 Prediction of stable delamination growth 
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Fig 10 R-curves for delamination growth in unidirectional double cantilever beam specimens 
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Fig 12 Schematic showing rate of delamination growth between ^tigue threshold 
and static interlaminar fracture toughness 
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Fig 13 Strain energy release rate analysis of instability related delamination growth 
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g 14 Cyclic compression fatigue life trends for T300/5208 graphite epoxy laminates 
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Fig 17 Edge delamination test for mixed-mode interlaminar fracture toughness 
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Fig 18 Interlaminar fracture toughness of graphite composites with different matrix resins 










